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We use statistical inference theory to explore the constraints from future galaxy weak lensing 
(cosmic shear) surveys combined with the current CMB constraints on cosmological parameters, 
focusing particularly on the running of the spectral index of the primordial scalar power spectrum, 
as. Recent papers have drawn attention to the possibility of measuring as by combining the 
CMB with galaxy clustering and/or the Lyman-a forest. Weak lensing combined with the CMB 
provides an alternative probe of the primordial power spectrum. We run a series of simulations with 
variable runnings and compare them to semi-analytic non-linear mappings to test their validity for 
our calculations. We find that a "Reference" cosmic shear survey with fsky = 0.01 and 6.6 x 10* 
galaxies per steradian can reduce the uncertainty on Us and by roughly a factor of 2 relative 
to the CMB alone. We investigate the effect of shear calibration biases on lensing by including 
the calibration factor as a parameter, and show that for our Reference Survey, the precision of 
cosmological parameter determination is only slightly degraded even if the amplitude calibration 
is uncertain by as much as 5%. We conclude that in the near future weak lensing surveys can 
supplement the CMB observations to constrain the primordial power spectrum. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 

The recent observations of the cosmic microwave 
background (CMB) from the Wilkinson Microwave 
Anisotropy Probe ( WMAP) mission confirmed the stan- 
dard cosmological model to a very high degree of accu- 
racy Q. The WMAP data alone is consistent with a 
spatially flat model dominated by dark energy and dark 
matter, with nearly scale-invariant, adiabatic and Gaus- 
sian initial perturbations consistent with the simplest in- 
flationary models |^ ^ i^J . 

One of the remaining outstanding issues is whether the 
shape of the primordial power spectrum is consistent with 
the theoretical predictions. This is usually expressed in 
terms of the slope of the primordial scalar spectrum 
and its running = dus/dlnk. Most of the existing 
models predict that the slope is close, but not identical, 
to scale invariant, rig 1, and that running is small, 
tts ~ (1 — n^)^. The best way to observationally settle 
this question is by combining data over a wide range of 
scales. Combining data from WMAP, the Cosmic Back- 
ground Imager (CBI) and the Arcminute Cosmology 
Bolometer Array Receiver (ACBAR) together with 
constraints from the Two Degree Field Galaxy Redshift 
Survey (2dFGRS) Ref. ^ found 0.05 error on and 
0.025 on as, with scale invariant model fitting the data 
at 1.3 a. Ref. 0] added Lyman-a forest constraints from 
[^13 and reduced the error on running to 0.016 (as well 
as found evidence for running at 1.9a), but Ref. 
showed that the assumed errors were too small and the 
current Lyman-a forest constraints do not add signifi- 
cantly to the constraints from the CMB. The constraints 
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on Us in [l| were obtained assuming no tensor modes, the 
strong degeneracy between running and tensors in the 
current data complicates the interpretation of the one- 
dimensional marginalized probability distribution for a^ 
113. Other authors [H Hi El [Hi who did not use 
Lyman-a forest, also find no significant evidence for a 
running spectral index, with errors at the level of 0.03. 

While the current data show no evidence for running, 
the errors are still large and will be improved in the fu- 
ture. The most immediate improvement will come from 
the Sloan Digial Sky Survey (SDSS),15]. The galaxy 
power spectrum analysis is not expected to improve the 
2dF results in a statistical sense (but will be an impor- 
tant check of systematics) , since on scales smaller than 
a few Mpc the nonlinear bias becomes intractable. On 
the other hand, the Lyman-a forest spectrum and bis- 
pectrum analysis of several thousand quasar spectra in 
SDSS will lead to an order of magnitude improvement of 
the Lyman-a constraints 16] . In combination with the 
CMB this is expected to give an error of 5 x 10^"^ on ag. 

There are however uncertanties associated with the 
baryonic physics of the Lyman-a forest. Future high-£ 
CMB data from the Planck satellite (and, to a lesser ex- 
tent, measurement of the second and third acoustic peaks 
from future WMAP data) will enable precision measure- 
ment of tts independently of the galaxy and Lyman-a 
data, but contamination from secondary anisotropics and 
foregrounds at high £ may be significant, since we are in- 
terested in effects of only a few percent. (The measure- 
ment of the first acoustic peak by WMAP is compara- 
tively clean due to the large cosmological signal.) 

It is therefore desirable to have yet another probe of the 
primordial power spectrum at high wavenumbers. Weak 
lensing is an attractive candidate since most of the poten- 
tial systematics are related to the difficulty of the obser- 
vations and are unrelated to the systematics associated 
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with the other cosmological probes we have discussed. 
The purpose of this paper is to explore the constraints 
from future galaxy weak lensing (cosmic shear) surveys 
combined with the current CMB results on cosmological 
parameters, particularly the running of the scalar spec- 
tral index. Weak lensing is one of the most promising 
tools for an era of precision cosmology as it probes di- 
rectly the distribution of the gravitating matter and does 
not suffer from the uncertainties associated with nonlin- 
ear bias or gas physics (for reviews, see |0, 0, and 
references therein). 

A few authors have used statistical inference theory to 
report forecasts on how well cosmic shear surveys com- 
bined with CMB data are able to constrain cosmological 
parameters. See, for example [23, Our goal is to 

extend these studies to include the running, as well as 
investigate the effects of possible systematics on the re- 
sults. More recently, Refs. combined CMB data 
from various experiments with the Red-Sequence Clus- 
ter Survey (RCS) weak lensing survey data (53 deg^) 
[23. [2^. Even with a modest survey, they found that 
their analysis reduces the uncertainties by comparable 
amounts or better than when the CMB results are com- 
bined with other types of probes because it breaks de- 
generacies present in the CMB data. Their results are 
consistent with the earlier forecasts. Other recent pa- 
pers that reported statistical forecasts on how well weak 
lensing will be able to constrain various cosmological pa- 
rameters include mHHHIlilSllllllllllll. 

The potential of weak lensing as a cosmological probe 
has been realized quickly and there are many ongoing, 
planned, and proposed surveys, such as the Deep Lens 
Survey {http://dls.bell-labs.com/) [H; the NO AO Deep 
Survey {http://www.noao.edu/noao/noaodeep/); the 
Canada- France-Hawaii Telescope (CFHT) Legacy Survey 
(http://www.cfht.hawan.edu/Science/CFHLS/) the 
Panoramic Survey Telescope and Rapid Response Sys- 
tem (http://pan-starrs.ifa.hawah.edu/); the Su pernova 
Acceleration Probe {SNAP; http://snap.l bl.gov/^ IHIll 
Ispl; and the Large Synoptic Survey Telescope (LSST; 
^ttp://www. lsst.org/lsstJiome. html) 40]. In this paper 
we will discuss parameter constraints possible with a Ref- 
erence Survey obtaining the shapes of n = 6.6 x 10® galax- 
ies per steradian over a fraction f^ky ~ 0.01 of the sky 
(413 deg^), with a peak redshift Zp = 1. We also explore 
the effect of varying fsky and n. The Reference Survey is 
somewhat more ambitious than the CFHT wide synoptic 
survey (172 deg^); roughly comparable to the SNAP wide 
survey; and significantly less ambitious than the LSST. 



II. FORMALISM AND MODEL 

A. Model parameters 

We consider the following basic parameter set for weak 
lensing: flmh^, the physical matter density; fl^, the frac- 
tion of the critical density in a cosmological constant; 



ns(fco) and a^, the spectral index and running of the pri- 
mordial scalar power spectrum at ko; Cg", the amplitude 
of linear fluctuations; Zp, the characteristic redshift of the 
source galaxies (see Eq. 0J; and and Cr, the calibration 
parameters (see ijll C|) . In order to combine this with in- 
formation from the CMB we include ilfc/i^, the physical 
baryon density; r, the optical depth to reionization; T/S 
the scalar-tensor fluctuation ratio. We assume a spa- 
tially flat Universe with flm + = 1, thereby fixing flm 
and Hq as functions of our basic parameters, and we do 
not include dynamical dark energy, massive neutrinos, 
or primordial isocurvature perturbations. The primor- 
dial power spectrum of scalar density fluctuation is given 
by M 
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We as- 
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— 0, i.e. we ignore higher-order terms 



in the Taylor expansion of the primordial power spec- 
trum. We have taken our pivot wavenumber to be 
fco — 0.05 h/Mpc. We use as fiducial model the best fit for 
the WMAP(lyr)-)-ACBAR-|-CBI data from the Markov 
Chains in Ref. Jii]: f^bft.^ = 0.0228, ri^/i^ = 0.139, 
r^A = 0.74, Us = 0.95, tts = (with ±0.04 variations), 
cTg = 0.85, T = 0.177, T/S = 0.265, Zp = 1.0, Q = 0.0, 
and C,r = 0.0. 



B. Convergence power spectrum 

We use the Limbe r ap proximation to the convergence 
power spectrum [42, |43, |4J| , which is valid on small an- 
gular scales ^ ^ 1 for which the gravitational lensing 
deflection can be approximated as a random walk due to 
many independent structures along the line of sight: 
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Here P^jj is the 3D nonlinear power spectrum of the 
matter density fluctuation, 6p/ p; x is the radial comov- 
ing coordinate; a{x) is the scale factor; and sinx X = 
K-^^^sm{K^^^X) is the comoving angular diameter dis- 
tance to X- The weighting function g{x) is the source- 
averaged distance ratio given by 
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^" I „ smg(x^ -x) , 

n{x ) ■■ TT^—dx , (3) 

smx X 



where n{x{z)) is the source redshift distribution normal- 
ized by j dz n{z) = 1. 

For our weak lensing calculations, we used the BBKS 
linear transfer function appropriate for cold dark 
matter with adiabatic initial perturbations. No baryonic 
correction was applied, since Vliyh^ is tightly constrained 
by the CMB and since the baryonic oscillations that ap- 
pear in the full transfer function are smoothed out by 
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projection effects wlien tlie convergence power spectrum 
is determined. We have used the analytic approximation 
to the growth factor of Ref. ji^. We used the recent 
nonhnear mapping procedure halofit to compute 
the non-hnear power spectrum. This procedure, as de- 
scribed in is based on a fusion of the halo model 
and an HKLM scaling |43| and is more accurate than 
the commonly used Peacock-Dodds mapping We 
discuss it further in section IIIII and check it versus nu- 
merical simulations. We contrast in Fig. ^ two curved 
convergence power spectra for as = —0.04, 0.04 and dis- 
play the sample variance errors averaged over bands in £ 
and the noise in the measurement. For n(z), we take the 
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FIG. 1: Curved convergence power spectra for the Reference 
Survey with Os — —0.04 (solid curve) and as = 0.04 (dashed); 
ko — 0.05h/Mpc and all other parameters are fixed at their 
fiducial values. The dotted line represent the noise power 
spectrum. The Reference Survey is sampling variance-limited 
out toix 2000. 



fitting function [SC 
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which peaks at Zp — 2zo. When we consider tomography, 
the sources are split in two bins: z < Zp and z > Zp. 
The normalized source redshift distributions for these two 
bins are: 
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for Zp < 2zo 
for Zp > 2z(). 
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systematically over- or under-estimated by a multiplica- 
tive factor. Physically, the principal source of this bias 
is incomplete correction for the circularization of images 
of galaxies by the point-spread function (PSF), although 
there are also noise- and selection-related contributions. 
The shear calibration bias is particularly dangerous be- 
cause the usual systematic error tests applied in weak 
lensing - e.g. decomposition of the shear field into E 
and B modes, cross-correlation of the shear map against 
maps of the PSF, etc. - are completely insensitive to 
this bias. Indeed, shear calibration bias mimics an over- 
all rescaling of the shear power spectrum, and proposals 
to circumvent it have thus far been based on detailed 
simulations of the observation. 

We have parameterized the shear calibration bias here 
using the power calibration parameter that is, the 
measured convergence power spectrum Pn{t} is given by: 



Pk(^) = (i + C.)PkW, 



(6) 



where Pk{£) is the power spectrum obtained in the ab- 
sence of calibration errors. Note that Cs refers to the 
calibration error of the power spectrum, which is twice 
the calibration error of the amplitude because the power 
spectrum is proportional to amplitude squared. 

When we consider tomography, we must also consider 
the relative calibration between the two redshift bins. 
This error affects the measured power spectrum [i) in 
accordance with: 
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where Ja ~ 0.32 is the fraction of the source galaxies in 
bin A and Jb ~ 0.68 is the fraction of the source galaxies 
in bin B. Mathematically, this means that if C,r — 0.01, 
then the power spectrum calibrations of the two redshift 
bins are offset 1% relative to each other, but the calibra- 
tion is correct when averaged over all the source galaxies. 
Physically, C,r parameterizes a redshift-dependent calibra- 
tion bias, which could arise e.g. from an incomplete PSF 
correction whose residual error depended on the galaxy's 
angular size or magnitude (which correlates with red- 
shift). 

It is possible for a redshift-dependent calibration bias 
to occur even without tomography, and hence we have 
included the relative calibration bias C,r even in our no- 
tomography parameter forecasts. In this case, the mea- 
sured convergence power spectrum is computed from: 



C. Calibration parameters 

Weak lensing surveys today are subject to various 
systematic errors that are comparable to the statisti- 
cal errors. One example is the shear calibration bias 
[siL Is^. fssL Is^. Isslls^ ■ in which the gravitational shear is 



P. w = fip. w + VAfsPr w + fsPr w , (8) 

where the P^^'' {€) are given by Eq. iQ. In principle the 
relative calibration bias can influence parameter estima- 
tion because it alters the effective source redshift distri- 
bution. However, we have found that it only slightly af- 
fects our parameter forecasts in the no-tomography case. 



4 



In a real experiment, (s and (r are not completely 
unknown, but rather are parameters of the experiment 
that can in principle be determined by simulating obser- 
vations. We have therefore imposed Gaussian priors of 
width CTpr(Cs) and crpr(Cr) on the calibration parameters 
(for simplicity we have taken C.s and C.r to be uncorrelated 
in the case of tomography); here the apr are the uncer- 
tainties in the shear calibration of the experiment. The 
prior curvature matrix (with diagonal elements cor- 
responding to the shear calibration parameters) is used 
in computing parameter uncertainties in accordance with 
Eq. (112). 

It is important to note that {Cs, t^r} is far from a com- 
plete parameterization of systematic errors. Other effects 
include spurious power from non-circular PSF, intrinsic 
alignments of galaxies, etc. We have not investigated 
these here, although clearly they must be minimized and 
the residual errors estimated if weak lensing is to evolve 
into a precision cosmological tool. 



III. NUMERICAL SIMULATIONS 

The HALOFIT code described in is a physically mo- 
tivated fitting formula to the results of numerical simu- 
lations. It is not obvious whether it is accurate enough 
for the present purpose: computing the derivatives of 
the non-linear power spectrum that are used to form 
the Fisher matrix. Ref. [i^l effectively varied all the 
parameters we use here, but their sampling of parame- 
ter space may not have been sufficiently dense to con- 
strain the interpolation near our point of interest. We 
performed a set of particle-mesh TV-body simulations to 
check HALOFIT in this context, using the Tree-Particle- 
Mesh (TPM) code described in [H^]- Each of our parame- 
ters is explicitly varied around the central model and the 
changes in the power spectra are compared to the predic- 
tions of the fitting formula (the figures we show are for the 
power at z = 0.5). The central model and variations are 
as = 0.85±0.05, fl^ = 0.74±0.04, fi^/i^ = 0.139±0.01, 
n = 0.95 ± 0.05, and a = 0.0 ± 0.04. We always compare 
AP/P, the fractional change in power, because most of 
the statistical error in AP is removed this way and the 
convergence with resolution is better. When we check 
how the simulations affect the Fisher matrix, we compute 
it using the formula AP = {APsim/ Psim)Pfit, where 
Psim means power from the simulation, and Pfn means 
power from halofit. 

Figure [21 shows the effect of changing a, for a pivot 
point fco = 0.2/i/Mpc that roughly corresponds to the 
scale of as and weak lensing. To be sure that our re- 
sults have numerically converged, we compare simula- 
tions with box size L — 320/i^^Mpc and N — 256"^ parti- 
cles (our mesh for the force calculation is always a factor 
of 2 finer than the mean particle spacing) to simulations 
with L = 160/i^^Mpc and N — 128^ (the same resolution 
as the bigger box), N = 256^, or = 512^. In Fig. Elwe 
note first that the difference between the two box sizes 




k [(Mpc/h)-"] 



FIG. 2: The fractional difference between the nonlinear mass 
power spectra ai z = 0.5 for a{ko — 0.2/i/Mpc) = 0.04 and 
a{ko = 0_.2ft/Mpc) = -0.04, [P{k,a = -f0.04) - P{k,a = 
-0.04)]/P(fc), where P{k) = [P{k,a = +0M) + P{k,a = 
— 0.04)]/2. All other parameters are fixed at their fiducial 
values. The dotted line is linear theory The solid line is from 
an L = 320^~^Mpc, A'^ = 256^ simulation. The open squares, 
filled squares, and triangles are from L = 160/i~^Mpc simula- 
tions with N = 512^, N = 256^, and N = 128^, respectively. 
Note that the triangle at the highest k is well beyond the 
resolution limit of the simulations. The dashed line is from 

HALOFIT. 



is probably insignificant (solid line vs. triangles). The 
L = 160/i~^Mpc, N = 256^ box has sufficient resolution 
to compute derivatives accurately out to fc = 10/i/Mpc 
(open vs. closed squares). Turning to the comparison 
between halofit and the simulations, we see that the 
fitting formula does well if one is not interested in the 
fine details of the nonlinear power; however, a Fisher 
matrix constraint that actually relies on the value of the 
derivative at fc ~ 0.2/i/Mpc (rather than just the fact 
that it is close to zero) would be inaccurate. We show in 
13 how the differences in the sc halofit and n-body mass 
power spectra are transferred to the convergence power 
spectra. 

For the combination of CMB and lensing we use the 
pivot point fco ~ 0.05/i/Mpc. Figure 0] shows the com- 
parison between the simulations and halofit for this fcg. 
The result is similar to the fco = 0.2/i/Mpc case, i.e., the 
basic trend is correct but the details are not completely 
correct. The same could be said of the results for varia- 
tion in n, shown in Fig. [S] 

The comparisons for variations in erg, and flmh^ 
are shown in Figures IHl and |S| respectively. Gener- 
ally the agreement for these cases seems better (i.e., very 
good). For fc ^ 2/i/Mpc the error in the prediction of the 
derivatives by the fitting formula is not large, even as a 
fraction of the values of the derivatives. 

Note that our comparison of halofit to the simu- 
lations has been restricted to the fractional change be- 
tween models, and implies nothing about the accuracy 
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FIG. 3: The fractional difference between the conver- 
gence power spectra for a{ko — 0.2/i/Mpc) = 0.04 and 
a{ko = 0^2/i/Mpc) = -0.04, [P4l,a = +0.04) - P4l,a = 
-0.04)]/P,(0, where P^l) = [P4l,a = +0.04) + P4l,a = 
— 0.04)]/2. All other parameters are fixed at their fiducial 
values. The solid line is from an L = 320/i"^Mpc, = 256^ 
simulation. The open squares and filled squares are from 
L = 160/i"^Mpc simulations with N = 512^ and N = 256^ 
respectively. The dashed line is from halofit. 
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FIG. 4: The fractional difference between the nonlinear mass 
power spectra at z = 0.5 for a{ko — O.OS/i/Mpc) = 0.04 and 
a{ko = 0.05/i/Mpc) = —0.04. The dotted hue is linear theory. 
The solid line is from an L = 320Mpc//i, N = 256^ simula- 
tion. The squares and triangles are from L = 160Mpc//i 
simulations with TV — 256'' and N = 128^, respectively. The 
dashed line is from halofit. 

of HALOFIT for predicting the absolute power spectrum, 
which is necessary when interpreting data. Our simula- 
tions are not suited to predicting the absolute power - 
on large scales they have big statistical fluctuations due 
to limited box size, while on small scales the power is 
suppressed by the limited PM resolution. Both of these 
effects cancel neatly in the fractional derivatives that we 
show, and for our Fisher matrix calculation this is suffi- 
cient. 
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FIG. 5: The fractional difference between the nonlinear mass 
power spectra sd, z — 0.5 for n — 1.0 and n — 0.9. The dotted 
line is linear theory. The solid line is from an L = 320Mpc//i, 
N — 256'' simulation. The squares and triangles are from 
L = laOMpc//!. simulations with iV = 256^ and N = 128^, 
respectively. The dashed line is from halofit. 
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FIG. 6: The fractional difference between the nonlinear mass 
power spectra a.t z = 0.5 for as = 0.9 and erg — 0.8. See Fig. 
|S|for line meanings. 

IV. PARAMETER FORECASTS 

A. Fisher-matrix analysis 

The uncertainty in the observed weak lensing spectrum 
is given by: 

where fsky = 8^7r/129600 is the fraction of the sky cov- 

1 /2 

ered by a survey of dimension 8 in degrees, (jf^t) ~ 
0.4 is the intrinsic ellipticity of galaxies. We assume for 
a Reference Survey a sky coverage of fsky = 0.01 and an 
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FIG. 7: The fractional difference between the nonhnear mass 
power spectra sd, z — 0.5 for f^A ~ 0.78 and SIa = 0.7. See 
Fig. |S]for hne meanings. 
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FIG. 8: The fractional difference between the nonlinear mass 
power spectra aX z — 0.5 for flmh^ = 0.149 and flmh? = 
0.129. See Fig. I^for line meanings. 



average galaxy number density of n « 6.6 x lO^sr"^; we 
will also investigate the effects of varying both of these 
parameters. 

The Fisher-matrix formalism for cosmological param- 
eter forecast has been proven in previous studies to be a 
powerful tool for estimating the statistical errors achiev- 
able by experiments [23, Ha, [sl, |63 . If the conver- 
gence field is Gaussian, and the noise is a combination 
of Gaussian shape and instrument noise with no intrinsic 
correlations, the Fisher matrix is given by: 



E 



{i+l/2)Uy dP^dP, 
(P. + hit) InY dv^ dv^ ■ 



(10) 



we have used £max — 3000 since on smaller scales, the as- 
sumption of a Gaussian shear field underlying Eq. H1U|) 
and the halofit approximation to the nonlinear power 



spectrum may not be valid. (Note that the shear field is a 
projection through many nearly independent structures, 
so by the Central Limit Theorem it can be well-described 
by Gaussian statistics even when the density perturba- 
tions are not Gaussian. Even with the Central Limit 
Theorem, deviations of the shear field from Gaussian 
statistics become large at the 1-2 arcmin scale IsHfS^.) 
On small scales, there is cosmological information in the 
small-scale non-Gaussianity (e.g. skewness) of the lens- 
ing field but we do not investigate this here. For the 
minimum £, we take the fundamental mode approxima- 
tion: 



360 deg 

e 



I sky 



(11) 



i.e. we consider only lensing modes for which at least 
one wavelength can fit inside the survey area. The sur- 
vey contains some information on larger angular scales, 
and for this reason the approximation Eq. (|ll|l might 
be considered conservative. However, it should also be 
noted that some of the planned surveys will scan discon- 
nected regions of the sky in order to provide additional 
systematic error checks, in which case the effective fmin 
is increased. 

The statistical error on a given parameter pa is then 
given by: 



a2(pJ«[(F-)-n)-i], 



(12) 



where H is the prior curvature matrix. We only impose 
priors on the source redshift and on the calibration pa- 
rameters (sec i)II C|) . For the Reference Survey, we take 
priors of cr(Cs) — criQr) — 0.02 on the calibration param- 
eters ^3 and (j{zp) — 0.05 on the source redshifts. Since 
the primary CMB anisotropies are generated at much 
larger comoving distance than the density fluctuations 
that give rise to weak lensing, it is a good approxima- 
tion to take them to be independent; in this case, we can 
add the Fisher matrices from lensing and CMB to yield 
combined constraints on cosmological parameters. This 
combination leads to significant improvements in param- 
eter estimation as we show in Fig. 



B. Adding tomography 



Tomography has been shown to improve significantly 
the measurements of cosmological parameters l64| . 
We explore the impact on parameter estimation from 
tomographic separation of the source galaxies into two 
redshift bins; see Eq. (Q. The Fisher formalism is gen- 
eralized here using: 

F^p = X:(^ + l/2)/.kyTr (^C7ipC7i|^) , (13) 



where is the covariance matrix of the multipole mo- 
ments of the observables " = P^'^ + N^'^ with N^'^ 
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FIG. 

with 
here. 



9: Improvements of the CMB parameter estimation 
weak lensing. The effect of increasing fsky is displayed 



the power spectrum of the noise in the measurement. 
These read 



pAA pAB 
pAB pBB ) ) 



and 
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(14) 



(15) 



The convergence power spectra for the two bins and their 
cross correlation are shown in Fig. 1101 
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FIG. 10: Convergence power spectra for bin A, bin B and 
their cross correlation, as — —0.04 and fco — 0.05h/ A4pc. 



and then combine these errors with the CMB analysis. 
We used the CMB parameter estimation uncertainties 
of Ref. ^3 (the correlation matrix was provided by 
the authors). As described in their paper, the authors 
have used a standard Monte Carlo Markov Chain us- 
ing the WMAP+CBI+ACBAR data (we will refer to 
these as "CMB"). In Table U we show the values for 
the uncertainties from CMB alone and then from a com- 
bination of CMB with weak lensing with and without 
tomography for our reference survey {fsky = 0.01 and 
n = 6.6 X lO^sr-i.) 

We summarize in Table Inland Table ITlll the results re- 
spectively for increasing values of fsky and n. We show 
the full correlation matrices for combined CMB and weak 
lensing observations with and without tomography in Ta- 
ble llVI In Table IVl we have fixed the calibration parame- 
ters and the characteristic redshift of the source galaxies. 
This significantly improves the precision of cosmological 
parameter estimates; we can see from the correlation ma- 
trix fTable lrV|) that most of this improvement comes from 
the fixing the source redshift. As can be seen from Table 
im very stringent cosmological constraints can in princi- 
ple be obtained for surveys covering thousands of square 
degrees; however, equally stringent controls over system- 
atic errors would be required. An additional caveat is the 
possibility of errors in the parameter constraints at very 
high fsky (i.e. the bottom row of Tables IITI and IV|l . since 
there the Fisher matrix is poorly conditioned and inaccu- 
racies in our computation of the derivatives 9Pk {£) / dp" 
are thus magnified. 

It is worth noting from the correlation matrix (Table 
IIV|) that the calibration parameters are not degenerate 
with any of the cosmological parameters considered. This 
lack of degeneracy is good news because it means that 
weak lensing surveys can be used for precision cosmology 
even if the absolute shear calibration cannot be deter- 
mined with high accuracy. Note that the current data 
analysis methods for lensing are estimated to have am- 
plitude calibration correct at roughly the ^ 5% level; the 
power calibration is twice this, or ±0.1 If we impose 
only this ±0.1 prior on (s and C,. (and retain the ±0.05 
prior on Zp), we see that the uncertainties on cosmologi- 
cal parameters are essentially unchanged (TablcQJ. With 
this weaker prior, the correlation of the calibration pa- 
rameters with Zp, CTg, and J^a becomes larger {\p\ > 0.5), 
and with even weaker priors the uncertainties in these 
parameters are degraded. 



DISCUSSION 



C. Results 

We generated the convergence power spectra using a 
weak lensing code that includes the formalism and fea- 
tures described in We derived parameter uncer- 
tainties for weak lensing with and without tomography 



In accord with previous studies, we find that when 
weak lensing surveys are combined with the CMB re- 
sults, the uncertainties on the cosmological parameters 
are reduced. The improvement can be up to an order of 
magnitude, depending on the size and depth of the sur- 
vey (see Table ^ and Table IIII|I . It is well known that 
weak lensing and CMB have different types of degenera- 
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TABLE I: A comparative table for paxameter estimation errors from our Reference Survey with fsky = 0.01 and n = 6.6 x 
10*sr"^ The first row (C) are uncertainties from the CMB= WMAP(lyr)+CBH-ACBAR alone. The following rows show 
combined errors for CMB plus weak Icnsing in the Reference Survey (CW), and CMB plus weak lonsing with tomography 
(CT). The lower-case suffixes a,b,c indicate the priors used: (a) prior of 0.1 on (s and (if applicable), and 0.05 on Zp; (b) 
prior of 0.02 on (a and fr) and 0.05 on Zp; and (c) prior of 10^'' on (Js, (^r, and Zp (i.e. effectively perfect knowledge of calibration 
and source redshift is assumed). For the rows labeled with [S], the A'^-body simulation was used to provide the matter power 
spectra instead of using halofit. Note that the simulated and halofit results are very similar. 







fT(n6/i^) 


<t{Qa) 


iT(cr8) 


(T{ns) 


o'(as) 


air) 


a{T/S) 








c 


0.013 


0.0012 


0.054 


0.083 


0.036 


0.039 


0.035 


0.154 








CWa 


0.0036 


0.0010 


0.0191 


0.024 


0.021 


0.017 


0.033 


0.143 


0.042 


0.092 


0.099 


CWb 


0.0030 


0.0010 


0.0167 


0.023 


0.021 


0.017 


0.033 


0.143 


0.041 


0.020 


0.020 


CWc 


0.0027 


0.0010 


0.0145 


0.022 


0.021 


0.017 


0.033 


0.142 


0.0001 


0.0001 


0.0001 


CTa 


0.0034 


0.0010 


0.0170 


0.021 


0.020 


0.016 


0.031 


0.140 


0.029 


0.075 


0.024 


CTb 


0.0029 


0.0009 


0.0159 


0.017 


0.019 


0.015 


0.030 


0.138 


0.027 


0.020 


0.015 


CTc 


0.0027 


0.0008 


0.0069 


0.010 


0.016 


0.013 


0.023 


0.130 


0.0001 


0.0001 


0.0001 


CWb[S] 


0.0032 


0.0010 


0.0169 


0.024 


0.022 


0.019 


0.034 


0.142 


0.045 


0.020 


0.020 


CWc[S] 


0.0027 


0.0010 


0.0151 


0.023 


0.022 


0.019 


0.034 


0.142 


0.0001 


0.0001 


0.0001 


CTb[S] 


0.0030 


0.0009 


0.0147 


0.017 


0.015 


0.014 


0.030 


0.125 


0.024 


0.020 


0.014 


CTc[S] 


0.0025 


0.0008 


0.0069 


0.010 


0.014 


0.013 


0.022 


0.123 


0.0001 


0.0001 


0.0001 



TABLE IL Parameter estimation errors for weak lensing combined with CMB: fgky is varied from 10 ^ to 1.0. and n = 
6.6 X 10* sr-^ We use priors of 0.02 on and Cr and 0.05 on Zp. Tomography is added in the second part of the table. 



fsky 






o-(nA) 


cricTa) 






cT(r) 


a{T/S) 


a{zp) 




<T{Cr) 


0.00001 


0.0052 


0.0011 


0.0266 


0.035 


0.034 


0.038 


0.035 


0.148 


0.050 


0.020 


0.020 


0.00010 


0.0037 


0.0011 


0.0216 


0.026 


0.032 


0.035 


0.035 


0.147 


0.050 


0.020 


0.020 


0.00100 


0.0033 


0.0010 


0.0182 


0.025 


0.025 


0.025 


0.034 


0.145 


0.048 


0.020 


0.020 


0.01000 


0.0030 


0.0010 


0.0167 


0.023 


0.021 


0.017 


0.033 


0.143 


0.041 


0.020 


0.020 


0.10000 


0.0028 


0.0009 


0.0160 


0.018 


0.019 


0.015 


0.031 


0.138 


0.026 


0.020 


0.020 


1.00000 


0.0022 


0.0008 


0.0132 


0.015 


0.016 


0.012 


0.026 


0.119 


0.018 


0.019 


0.020 


0.00001 


0.0051 


0.0011 


0.0264 


0.034 


0.034 


0.038 


0.035 


0.148 


0.049 


0.020 


0.020 


0.00010 


0.0037 


0.0011 


0.0215 


0.025 


0.031 


0.035 


0.034 


0.147 


0.047 


0.020 


0.020 


0.00100 


0.0032 


0.0010 


0.0180 


0.020 


0.024 


0.023 


0.032 


0.143 


0.037 


0.020 


0.019 


0.01000 


0.0029 


0.0009 


0.0159 


0.017 


0.019 


0.015 


0.030 


0.138 


0.027 


0.020 


0.015 


0.10000 


0.0026 


0.0008 


0.0122 


0.014 


0.016 


0.012 


0.026 


0.122 


0.019 


0.018 


0.008 


1.00000 


0.0019 


0.0006 


0.0067 


0.009 


0.010 


0.007 


0.020 


0.091 


0.011 


0.014 


0.003 



TABLE III: Parameter estimation errors for weak lensing combined with CMB: n is varied from n = 5.0 x 10®sr~^ to n = 
1.0 X lO^sr"^; fsky = 0.01. With a prior of 0.02 on i^^ and and 0.05 on Zp. Tomography is added in the second part of the 
table. 



n 






CT(nA) 




a{ns) 


(T{as) 


air) 


a{T/S) 


(T{Zp) 


<^{Cs) 


<T{Cr) 


5.0 X 10" 


0.0045 


0.0011 


0.0239 


0.032 


0.034 


0.038 


0.035 


0.149 


0.050 


0.020 


0.020 


1.0 X 10^ 


0.0035 


0.0011 


0.0194 


0.028 


0.030 


0.033 


0.035 


0.147 


0.049 


0.020 


0.020 


5.0 X 10^ 


0.0032 


0.0010 


0.0173 


0.024 


0.023 


0.022 


0.034 


0.145 


0.046 


0.020 


0.020 


1.0 X 10** 


0.0032 


0.0010 


0.0170 


0.024 


0.022 


0.020 


0.034 


0.144 


0.044 


0.020 


0.020 


5.0 X 10* 


0.0031 


0.0010 


0.0167 


0.023 


0.021 


0.017 


0.033 


0.143 


0.041 


0.020 


0.020 


1.0 X 10^ 


0.0030 


0.0010 


0.0166 


0.022 


0.021 


0.017 


0.033 


0.143 


0.040 


0.020 


0.020 


5.0 X 10" 


0.0048 


0.0011 


0.0247 


0.033 


0.034 


0.038 


0.035 


0.149 


0.050 


0.020 


0.020 


1.0 X 10^ 


0.0036 


0.0011 


0.0197 


0.028 


0.030 


0.034 


0.035 


0.147 


0.049 


0.020 


0.020 


5.0 X 10^ 


0.0032 


0.0010 


0.0173 


0.021 


0.022 


0.021 


0.033 


0.143 


0.040 


0.020 


0.019 


1.0 X 10* 


0.0031 


0.0010 


0.0168 


0.019 


0.021 


0.019 


0.032 


0.142 


0.035 


0.020 


0.018 


5.0 X 10* 


0.0029 


0.0009 


0.0159 


0.017 


0.019 


0.016 


0.031 


0.138 


0.027 


0.020 


0.015 


1.0 X 10^ 


0.0028 


0.0009 


0.0148 


0.015 


0.018 


0.014 


0.030 


0.133 


0.023 


0.019 


0.009 
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TABLE IV: Correlation matrix for parameter estimation errors for weak lensing combined with CMB. n = 6.6e8 and fsky = 0.01. 
With a prior of 0.02 on (r and Cs, and 0.05 on Zp. The above-diagonal elements correspond to a survey without tomography; 
the below-diagonal elements correspond to a survey with tomography. 











ns 


as 


Zp 


(s 






T 


T/S 




1 


-0.325 


-0.020 


0.135 


-0.266 


-0.430 


-0.137 


0.017 


0.378 


-0.533 


0.313 


nA 


-0.375 


1 


0.754 


0.667 


-0.583 


0.464 


0.122 


-0.018 


0.558 


0.728 


0.393 




-0.230 


0.877 


1 


0.632 


-0.679 


-0.160 


-0.091 


0.008 


0.541 


0.757 


0.375 


Us 


0.072 


0.651 


0.620 


1 


-0.774 


0.104 


0.022 


-0.002 


0.788 


0.393 


0.794 


as 


-0.228 


-0.550 


-0.626 


-0.770 


1 


-0.076 


0.026 


0.005 


-0.716 


-0.433 


-0.752 


Zp 


-0.386 


0.818 


0.570 


0.444 


-0.455 


1 


-0.052 


0.020 


0.092 


0.126 


0.096 


Cs 


-0.156 


0.078 


-0.217 


-0.028 


0.106 


-0.045 


1 


0.002 


0.014 


0.004 


0.013 


Cr 


0.086 


0.158 


0.386 


0.221 


-0.204 


-0.243 


0.083 


1 


-0.003 


-0.005 


-0.004 




0.349 


0.527 


0.516 


0.762 


-0.698 


0.371 


-0.033 


0.185 


1 


0.372 


0.673 


T 


-0.680 


0.726 


0.732 


0.329 


-0.321 


0.572 


-0.049 


0.176 


0.311 


1 


0.136 


T/S 


0.294 


0.355 


0.357 


0.786 


-0.756 


0.266 


-0.028 


0.133 


0.650 


0.073 


1 



TABLE V: Parameter estimation errors for weak lensing combined with CMB: Here we have fixed ^r, Cs ^^'i by imposing 
priors of width 10~* on these parameters; fsky is varied from 10~^ to 1.0; n = 6.6 x lO^sr"^. Tomography is added in second 
part of the table. 



fsky 






^(Oa) 


o-(cr8) 


cr(ns) 


cr(as) 


a(r) 


a{T/S) 




^(Cs) 




0.00001 


0.0050 


0.0011 


0.0260 


0.034 


0.034 


0.038 


0.035 


0.148 


0.0001 


0.0001 


0.0001 


0.00010 


0.0034 


0.0011 


0.0206 


0.025 


0.032 


0.035 


0.035 


0.147 


0.0001 


0.0001 


0.0001 


0.00100 


0.0029 


0.0010 


0.0169 


0.023 


0.025 


0.025 


0.034 


0.145 


0.0001 


0.0001 


0.0001 


0.01000 


0.0027 


0.0010 


0.0145 


0.022 


0.021 


0.017 


0.033 


0.142 


0.0001 


0.0001 


0.0001 


0.10000 


0.0026 


0.0009 


0.0105 


0.016 


0.017 


0.013 


0.027 


0.132 


0.0001 


0.0001 


0.0001 


1.00000 


0.0022 


0.0007 


0.0044 


0.007 


0.011 


0.008 


0.019 


0.104 


0.0001 


0.0001 


0.0001 


0.00001 


0.0049 


0.0011 


0.0258 


0.033 


0.034 


0.038 


0.035 


0.148 


0.0001 


0.0001 


0.0001 


0.00010 


0.0033 


0.0011 


0.0204 


0.024 


0.031 


0.035 


0.034 


0.147 


0.0001 


0.0001 


0.0001 


0.00100 


0.0029 


0.0010 


0.0146 


0.019 


0.023 


0.022 


0.031 


0.140 


0.0001 


0.0001 


0.0001 


0.01000 


0.0027 


0.0008 


0.0069 


0.010 


0.016 


0.013 


0.023 


0.130 


0.0001 


0.0001 


0.0001 


0.10000 


0.0024 


0.0007 


0.0024 


0.003 


0.012 


0.008 


0.019 


0.112 


0.0001 


0.0001 


0.0001 


1.00000 


0.0019 


0.0006 


0.0008 


0.001 


0.007 


0.004 


0.017 


0.085 


0.0001 


0.0001 


0.0001 
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cies in their parameters which are nicely broken when 
combined together. In particular, weak lensing does not 
suffer from the well-known angular diameter distance de- 
generacy; and it probes smaller comoving scales than the 
CMB, which means that it has a different degeneracy di- 
rection in the (us, as) plane. The improvements from in- 
cluding weak lensing are especially notable for as, ^mh^ 
and flA- 

Motivated by recent discussion concerning the running 
of the primordial power spectrum a^, we include it as a 
parameter in the weak lensing analysis. We find that 
for small surveys, modest improvement is obtained for 
Us and as- Our Reference Survey can reduce a{ns) and 
a{as) by roughly a factor of two. For surveys within the 
near future (see Refs. in SjJ), weak lensing can be used 
to provide complementary constraints to detect a possi- 
ble running spectral index, but may not be able to verify 
the results obtained by combined CMB-t-Lyman-a forest 
analysis, which should give another factor of 2-4 lower 
a (as). A detailed comparative study of the constraints 
from the two probes is left for future work. We find that 
tomography improves in particular the uncertainty on erg 
and Q\. For the Reference Experiment, we do not find 
large improvements in the other parameters from tomog- 
raphy except for the most ambitious surveys. This can be 
attributed to the strong parameter degeneracies present 
even when tomography is used (see Table llV|l . most no- 
tably the degeneracy between SIa and Zp. If the source 
redshift distribution is known accurately from a spectro- 
scopic redshift survey, then this degeneracy is lifted, and 
tomography becomes a powerful tool for measuring cos- 
mological parameters, including and (see Ta- 

ble 0. A comparison of Tables Hll and Ivl shows that the 
parameter constraints from tomography are more sensi- 
tive to the systematics parameters than the constraints 
without tomography, hence the benefits of tomography 
can only be fully realized if systematic errors are tightly 
controlled. 

We expanded the usual weak lensing parameter space 
to include two calibration parameters in addition to 
the characteristic redshift of source galaxies. We hope 
that the present analysis will encourage weak lensing ob- 
servers to expand their likelihood analyses to include a 
parameterization of systematic errors. It is already com- 
mon practice to report results marginalized over the char- 
acteristic source redshift j0 |3, US or to treat it as a 
systematic to be added in quadrature to statistical errors 
|65|. Ultimately it would be desirable to include not just 
the calibration factors and characteristic source redshift 
but also the full redshift distribution of sources, intrinsic 
alignments, etc. Additional data (such as spectroscopic 
redshifts) and detailed analysis and simulations may be 
required in order to constrain some of these parameters. 



However, in some cases it is possible to obtain informa- 
tion about the systematics parameters from the data it- 
self. For example, the nonlinear portion of the conver- 
gence power spectrum provides joint constraints on the 
shear calibration biases and Cr and the cosmological 
parameters. 

We have shown, at least at the level of our Refer- 
ence Survey, that halofit provides a sufhciently good 
fit to A^-body simulations for use in parameter forecast- 
ing studies. The largest discrepancy is 25% in the error 
on Us when tomography is included, but most of the er- 
ror estimates agree to better than 10%. Note that this 
does not imply that the nonlinear mapping is sufficiently 
accurate for analysis of Reference Survey data, since the 
best-fit values of cosmological parameters can be signifi- 
cantly affected by small errors in theoretical predictions 
even if the fractional change in the Fisher matrix is small. 

In summary, we have shown that weak lensing, supple- 
mented with the one-year WMAP data, has the poten- 
tial to very precisely measure the running of the scalar 
spectral index. This would provide a third measure- 
ment of as in addition to those using the Lyman-a forest 
and/or galaxy power spectrum, and the precision mea- 
surement of the high-^ CMB power spectrum expected 
from Planck. Weak lensing observations are currently 
progressing rapidly and the data required to significantly 
improve WMAP constraints on as should be available in 
the foreseeable future. By providing an additional and 
mostly independent measurement of the value of as (and 
more generally the scalar power spectrum) obtained by 
these other methods, weak lensing will help us to under- 
stand the spectrum of primordial scalar fluctuations in 
the universe and thus provide valuable information on 
the mechanism of their generation. 
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